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ABSTRACT 

We discuss a new scenario for the formation of intermediate mass black holes in dense star 
clusters. In this scenario, intermediate mass black holes are formed as a result of dynamical in¬ 
teractions of hard binaries containing a stellar mass black hole, with other stars and binaries. 
We discuss the necessary conditions to initiate the process of intermediate mass black hole 
formation and the influence of an intermediate mass black hole on the host global globular 
cluster properties. We discuss two scenarios for intermediate mass black hole formation. The 
SLOW and FAST scenarios. They occur later or earlier in the cluster evolution and require 
smaller or extremely large central densities, respectively. In our simulations, the formation of 
intermediate mass black holes is highly stochastic. In general, higher formation probabilities 
follow from larger cluster concentrations (i.e. central densities). We further discuss possible 
observational signatures of the presence of intermediate mass black holes in globular clusters 
that follow from our simulations. These include the spatial and kinematic structure of the host 
cluster, possible radio. X-ray and gravitational wave emissions due to dynamical collisions 
or mass-transfer and the creation of hypervelocity main sequence escapers during strong dy¬ 
namical interactions between binaries and an intermediate mass black hole. All simulations 
discussed in this paper were performed with the MOCCA Monte Carlo code. MOCCA ac¬ 
curately follows most of the important physical processes that occur during the dynamical 
evolution of star clusters but, as with other dynamical codes, it approximates the dissipative 
processes connected with stellar collisions and binary mergers. 
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1 INTRODUCTION 

Recent high resolution observations of globular clusters (GCs) have 
provided a detailed picture of their physical status, revealing com¬ 
plex phenomena connected with multiple stellar populations, bi¬ 
nary evolution and the Galactic tidal field. Despite such great ob¬ 
servational progress, many theoretical uncertainties remain in un¬ 
derstanding GC formation. These include the origins of their mul¬ 
tiple stellar populations, the properties of their primordial binary 
populations and the possible presence of both stellar mass and in¬ 
termediate mass black holes (IMBHs). 

IMBHs provide a missing link between stellar-mass black 
holes (BH) and supermassive BHs (SMBH). There are several def¬ 
initions of IMBHs available in the literature based on their charac¬ 
teristic mass range but, for the purpose of this paper, only IMBHs 
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with masses between about 100 — 10 5 Mq are considered. Obser¬ 
vational confirmation of the existence of IMBHs would have an 
important impact on a number of open astrophysical problems re¬ 
lated to the formation of SMBHs and their host galaxies, the ori¬ 
gins of ultraluminous X-ray sources in nearby galaxies and the 
detection of gravitational waves. In the Milky Way galaxy, the 
presence of IMBHs in the cores of some GCs has been debated 
for a long time. There are many theoretical and observational ar¬ 
guments in favou r of the formation of IM BHs in the centres of 
Galactic GCs fe.g. lLtitzgendorf et al.ll2013l . and reference therein), 
but no firm observational confirmation of their presence has been 
presented to date. Notwithstanding, several claims of their pres¬ 
ence have been made, based on applying Je ans’ model to detai led 
photometric and spectroscopic data ( e.g. iGebhardt, Rich & Hoi 
I2005L Ivan der Marel & Andersonl 1201 Ol . and references therein). 
More concretely, several direct observations of the X-ray emis¬ 
sion from distant star clusters have revealed ultraluminous X- 
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ray sources at off-ce ntre positions within their host galaxies (e.g. 
iMaccarone & Servillatl2008h . These observations suggest the pres¬ 
ence of IMBHs with masses of about a few thousand Mq. 

Several widely different theoretical scenarios have been pro¬ 
posed in the literature for the formation of IMBHs in GCs: (1) the 
formation of very massive P opulation II I stars and their subsequent 
collapse into a massive BH l iMadau & Reesl200lh . (2) the runaway 
merg ing of main sequence stars in young and very dense star clus - 
ters JPortegies Zwart et aPl2004l : iGiirkan, Freitag & Rasidl20o3) . 
(3) the accretion of residual gas o nto stellar mass BHs formed 
from first generation stars (Leigh et al.ll2013ah . A strong motiva¬ 
tion for searching for IMBHs is the observed relation between the 
blac k hole (BH) mass and the velocity dispersion of its host galaxy 
(e.g. Ferrarese,& Merritt 200Cil ;|Gtiitekin e t al]2009l . and references 
therein). Extrapolating this relation to the velocity dispersions char¬ 
acteristic of Galactic GCs, one expects to find central BHs with 
masses similar to those characteristic of IMBHs. 

Scenario (1) invokes the direct collapse of Population III stars 
formed out of unmagnetized metal-free gas, and implies that the 
initial mass function (IMF) was extremely top heavy. Stars with 
masses larger than 200 — 300Mq might have formed in so-called 
minihalos with masses about 10 5 Mq at redshifts z about 20. Af¬ 
ter only a few Myr, such extremely massive stars would have 
formed BHs with masses about 200 A/q , loosing only a small frac¬ 
ti on o f their mass in the process (Fryer. Wooslev & Hegetil200ll : 
I Whalen & Frvei|[20 1 2l) . These BHs could then have acted as the 
seeds for the formation of IMBHs in GCs. But, doubts remain re¬ 
garding the validity of this scenario. For instance, how could such 
a massive BH formed from a Population III star end up in a GC 
consisting of Population II stars? Recentlv. lBelczvnski et alj|20l4) 
discussed extensively the possibility of formation of very massive 
stars (with masses greater than 150Mq) in star-forming regions. 
They argue that such massive stars are currently forming, and that 
they could be IMBH seeds. 

Scenario (2) invokes the runaway merging of massive main se¬ 
quence (MS) stars. In order for this to occur, the time-scale for mass 
segregation of the most massive stars (about 100 Mq) has to be 
shorter than the stellar evolution time-scale for those stars. Other¬ 
wise, they will lose a substantial amount of mass before segregating 
to the cluster centre, where further mass-loss powers cluster expan¬ 
sion. Then, the core will be of low density and consist of relatively 
low-mass evolved stars, instead of being dense and populated by 
massive stars. Additionally, the velocity dispersion in a collapsed 
cluster cannot be larger than a few hundred km/s, otherwise col¬ 
lisions between massive MS stars will be disruptive and result in 
significant mass-loss. Regardless, if the right conditions are met, a 
very massive MS star can form from direct physical collisions be¬ 
tween many lower mass MS stars. The resulting collision product 
will shortly end its evolutio n as a very massive BH, an IMBH seed. 

Recently ILeigh et aD d2013al) proposed a new model for gas 
expulsion in very massive star clusters, which might relate to the 
formation of multiple stellar populations in GCs. The basic idea 
is connected with gas accretion onto BHs, formed from the most 
massive stars. The interaction between the BHs and the interstellar 
mediu m should spe edup and even enhance their mass segregation 
( ILeigh et al.ll2014a) . possibly increasing their masses in the pro¬ 
cess. This could lead to the rapid formation of a dense and (possi¬ 
bly) massive subsystem of BHs. Further evolution can lead to the 
formation of an IMBH, and/or significantly deplete GCs of their 
primordial gas reservoirs (via direct accretion onto BHs and/or me¬ 
chanical/radiative feedback). 

BHs cannot be observed directly. Indirectly, however, they are 


observable via radio. X-ray and gravitational radiation due to var¬ 
ious types of accretion processes. Observations of the radio and 
X-ray emissions from distant galaxies have led to detections of 
ultraluminous X-ray sources at off-centre positions, which possi- 
bly_are young star clusters, (e.g. Matsum oto et~a]l200ll : ISoria et aD 
1 20IQ . and references therein). These observations are indicative 
of the presence of IMBHs with masses of about a few thousand 
Mq. These masses are in agreement with what is expected from 
extend ing the observe d M - a correlation for supermassive BHs 
dMagorian et al.lfl998l) to lower galaxy masses. 

IMBHs may also be indirectly observable via subtle changes 
in the cluster structure in their vicinity. Theoretical investiga¬ 
tions and numerical simulations point out that if an IMBH is 
present in the centre of a GC the velocity dispersion profile 
(VDP) should rise in the centre, and the central surface brightness 


profile (SBP) should be shallow (e.g. Baumgardt. Makino & Hut 

2005 

: Trend, Vesperini & Pasquato 20ld: Novola & Baumgardt 

2011 

, and references therein). Thus, when hunting for IMBHs 


in GCs, the best candidates are massive GCs with shallow cen¬ 
tral SBPs and centrally rising VDPs. However, this straightfor- 
wa rd prediction sho uld b e ap plied with caution. As was shown 
by iBaumgardt. Makino & Hutl (2005j), the massive evolved stellar 
population in the core of a GC can also show a rising VDP. Con¬ 
versely, the presence of a massive stellar mass BH binary in the 
core can produce a shallow central SBP dHurlevll2007l ). 

In Section[2] we briefly describe the Monte Carlo code, called 
MOCCA, used in this paper to perform simulations of GC evolu¬ 
tion. In Section [3] we describe the initial conditions for our GC 
models. Next (Section [4]l we describe a new scenario for IMBH 
formation, and discuss the conditions needed for this to occur. The 
observational implications of our scenario for IMBH formation in 
GCs are discussed in Section[5] In the final two Sections, we sum¬ 
marise our results and discuss some limitations and future develop¬ 
ments of the MOCCA code. 


2 METHOD 

The MOCCA (MOnte Carlo Cluster simulAtor) code used for the 
star cluster simulations presented here is the Monte Carlo code 
based o n Henon’s implementation of the Monte Carlo method 
dllenonl Il97ll ). which was further substantially develop ed by 
Stodolkiewicz in the early eighties (StodolkiewiczlIl986l) . This 
method can be regarded as a statistical way of solving the Fokker- 
Planck equation. The basic assumptions behind the Monte Carlo 
method are: (1) spherical symmetry, which makes it easy to quickly 
compute the gravitational potential and stellar orbits at any place in 
the system. This is a very severe assumption, e.g. the evolution of 
a system with rotation cannot be investigated. In specified spheri¬ 
cally symmetric potentials every star is characterized by its mass, 
energy and angular momentum. Each star moves on a rosette or¬ 
bit. There is no need for integration of an orbit in the Monte Carlo 
method, since it is easy to calculate the position of each star along 
its orbit. In this way, the Monte Carlo method is very fast; (2) clus¬ 
ter evolution is driven by two-body relaxation, and the time step at 
each position in the system is proportional to the local relaxation 
time. This is the second reason why the Monte Carlo method is so 
fast. Generally, Monte Carlo codes cannot follow any physical pro¬ 
cesses with characteristic time scales much shorter than the local 
relaxation time. 

The MOCCA code jCiersz et al.11201 3l) is a heterogeneous 
code, composed of independent modules (e.g. single and binary 
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star evolution, small N-body scattering, etc.) that together model 
the entire cluster evolution. It is able to follow most of the impor¬ 
tant physical processes that occur during the dynamical evolution of 
star clusters. In addition to the inclusion of relaxation, which drives 
the dynamical evolution of the system, MOCCA includes: (1) syn¬ 
thetic stellar evolution of single ( SSE) and binary (BSE) star s us- 
ing the prescriptions provid ed bv lHurlev, Pols. & Toutl d200Ch and 
iHurlev. Tout, & Polsl 1200 2l) (from now on we will refer to both 
codes together as BSE), (2) direct integration pro cedures for small 
N sub-systems using the Fewbody code dFregeau et al 120041) . and 
(3) a realistic treatment of escape p r ocesse s in tidally limited clus¬ 
ters based on lFukushige & Heggiel d200Ct) . Here, the escape of an 
object from the system is not instantaneous, but delayed in time. 

The MOCCA code has b een extensively tes t ed against the 
results of N-bod y simulations dGiersz et alj 1201 3l : iHeggiel 120141 : 
[Long et alj|2015l . and references therein). The agreement between 
these two different types of simulations is excellent. This includes 
the global cluster evolution, mass segregation time scales, the treat¬ 
ment of primordial binaries (energy, mass and spatial distributions), 
and the numbers of retained neutron stars (NS) and BHs. The 
key assumption implemented in the MOCCA code that guaran¬ 
tees this agreement is that, throughout the entire cluster, the time 
scale for significant evolution is always a small fraction (about 
0.001 - 0.01) of the local relaxation time, but larger than the lo¬ 
cal crossing time. It is not possible to directly compare results 
from MOCCA and N-body for extremely large mass ratios. The 
only analogous N-body simulations presented in the literature are 
for small N and usually do not contain all the necessary physics 
- e.g . stella r/binary evolution, an artificially truncated IMF (e.g. 
iKonstantinidis, Am aro-Seoane. & Kokkotasl d2013l) ). However, the 
recent paper bv lLeigh et al. d2014bl) clearly shows that the results 
for binary star populations evolving in the presence of an IMBH 
are the same for N-body and MOCCA. This gives indirect evidence 
and some confidence that the chosen time-step in the MOCCA sim¬ 
ulations is appropriate. Also, to minimise density fluctuations asso¬ 
ciated with the most massive stars due to their small numbers, a 
special procedure was introduced to calculate the probability of bi¬ 
nary formation using a sm oother average density in the vicinity of 
the binaries dGierszIpOOll) . The relaxation of very massive stars via 
long-range encounters with other mostly low-mass stars, which is 
pivotal to IMBH formation in our simulations, is treated accord¬ 
ing to the standard two-body prescription given in lHeno a fi97ih . 
Since we do not directly calculate the relative force between every 
pair of particles at each time step, one might naively expect this 
approach to introduce some error into our results for IMBH for¬ 
mation. However, as is shown later in the paper (see Section [5T1 >. 
the spatial and kinematic structures of our model star clusters with 
IMBH formation reasonably reproduce the results of analogous N- 
body simulations performed by other authors. 

We would like to stress that the Fewbody code models only 
direct gravitational interactions between stars and binaries. The hy¬ 
drodynamics, gravitational radiation (GR), tidal effects (distortions 
and dissipation), mass-loss, etc. are not accounted for. Using the 
sticky-star approximation, we make the general assumption that, 
if the stellar radii overlap during an encounter, then the hydrody¬ 
namics will be sufficiently non-negligible to bind the stars at a 
sufficiently small separation that they are well within the stellar 
envelope of the collision product, which then subsequently settles 
back on to a standard stellar evolution track on a Kelvin-Helmholtz 
timescale (which should typically be longer than the total encounter 
duration, but shorter than the time until a subsequent encounter oc¬ 
curs), with the original components ultimately merging inside the 


envelope (rapidly - i.e. on a timescale much shorter than the en¬ 
counter time) to form the core of the new star. Thus, to first-order, 
we assume that after a collision occurs, the product is puffy for the 
remainder of the encounter, but returns to having a radius normal 
for a star of that mass post-encounter (if it avoids a subsequent col¬ 
lision). The subsequent stellar evolution of such a thermally relaxed 
star is treated by the BSE code. These assumptions are at least con¬ 
sistent with our current understanding of the relevant physics, and 
to address them further is beyond the scope of the present analy¬ 
sis. Notwithstanding, they should eventually be re-visited in more 
detail. 

Importantly, all of the neglected physics of significance to our 
results act as sources of energy damping during dynamical interac¬ 
tions. This includes tides and gravitational radiation. In general, we 
expect damping to facilitate collisions, and increase the frequency 
of collisions relative to our calculations. Similarly, in the case of 
stable or long-lived triples (which we neglect in our simulations), 
the presence of a distant third companion could reduce the time- 
scale for a subsequent encounter, increase the collis ion probability 
during individual encounters jLeigh & GelleJl2012h or even drive 
the inner binary t o shorter periods or even merge due to Kozai- 
Lidov oscillations ( 'Perets & Fabrvckvi2009b . Hence, our calculated 
collision rates should increase if these damping sources and stable 
triples are included in our analysis. Naively, we expect this to in¬ 
crease the rate of IMBH mass build-up, which will be the subject 
of future work. 

The accuracy and stability of the Fewbody code was exten¬ 
sively tested for very large mass ratios. For the mass ratios of inter¬ 
est to this paper, we have verified that the Fewbody code is able to 
successfully repr oduce the energy an d angu lar momenta exchange 
rates presented in lSesana, Haardt & Madaul d2015h . 

The MOCCA code provides as many details as N-body codes, 
while being much faster. It can follow the time evolution and move¬ 
ment of particular objects, and needs only about a day to complete 
the evolution of a realistic globular cluster. So, instead of just one 
N-body model, several hundreds of MOCCA models can be com¬ 
puted (with different initial conditions) in a comparable amount of 
time. Thus, the MOCCA code is ideal not only for dynamical mod¬ 
els of particular clusters, but also large surveys. 


3 INITIAL CONDITIONS 

The new scenario for IMBH formation presented in this paper ap¬ 
peared as a byproduct fr om two earlier projects. The aim of the 
first project l lLeigh et alj|2013bl) was to explain the observed cor¬ 
relations: between cluster mass and binary fraction - the lar ger the 
cluster mass, the smaller the binary fraction ( IMilone et al .1120120 : 
and mass function (MF) power-law index and cluster concentration 
- the more co ncentrated the cl uster , the larger the MF power-law 
index foe Marc hi. Paresce & Pulonell 20071) . The aim of the second 
nroiecti lLeigh et alj2015al) was to constrain the initial properties of 
primo rdial binaries by comparison with observations ( IMilone et al.l 
120121) and to check if star cluster simulati ons with initial conditions 
drawn from the invariant iKroupal i j 19951 ) distributions are able to 
recover the observed spatial distributions of binaries. All together, 
for these projects about 400 models of star clusters (SC) were sim¬ 
ulated by the MOCCA code. 

The parameters of models run for the aforementioned projects 
are listed in Table Q] Generally, more massive models have larger 
concentrations (measured as the ratio between tidal and half-mass 
radii C = Rt/Rh - i.e. not the standard concentration parame- 




























































4 M. Giersz, N. Leigh, A. Hypki, N. Liitzgendorf & A. Askar 


Table 1. Initial conditions for the MOCCA simulations. 


N 

Rt (pc) 

Rt/Rh 

fb 

CLmax (AU) 

IMF 

Kicks (km/s) 

N = 5.0xl0 4 
King Wo = 6 

38 

20, 25, 30 

50, 75 

0.1 

100 

Period 

IMF2, IMF3 

0, 265 
Fallback 

N = l.OzlO 5 

King Wo = 6 

69 

10, 20 

0.1 

100 

Period 

IMF3 

IMF2-m 

265 

Fallback 

N = 2.0xl0 5 

King Wo = 6 

69 

10, 20 

0.1 

100 

Period 

IMF3 

IMF2-m 

265 

Fallback 

N = 3.0a;10 5 

King Wo = 6 

35,50 

69 

35, 50, 60 

75, 100 

0.1, 0.3 
0.7, 0.95 

100,200 

400, Period 

IMF2, IMF3 

IMF2-m 

0, 265 
Fallback 

N = 1.8x10 h 

King Wo = 6 

125 

60, 75, 90 
100, 125 

0.1 

100 

Period 

IMF2, IMF3 

IMF2-m 

0, 265 
Fallback 

N = 5.5xl0 b 

N = 6.9a:10 5 
King Wo = 6 

181,98, 68 

50. 100 

0.95 

Period 

IMF2, IMF3 

265 

N = l.lxl0 H 

King Wo = 6 

220, 120, 75 

50, 100, 200 

0.95 

Period 

IMF2, IMF3 

265 

N = 6 . 62 IO 4 
King Wo = 6 

3.6 

6.8 

0.0, 0.1 
0.95 

100 

Period 

IMF2 

190 

Fallback 


Notes: BH and NS kicks are the same, except the case of mass fallback (iBelczynski, Kalogera & Bulik 


fraction, a ma x - maximum value for semi-major axis (distribution uniform in log(a)). 


2002 ). 


1993). 


1995). 


IMF2 - two segmented IMF jKroupaluOOll) , IMF3 - three segmented j Kroupa, Tout & Gilmore 
IMF2-m - modified by changing oq ow and othigh » Period - binary period distribution iKroup" 
fb - binary fraction, N - number of stars, Rt - tidal radius, - half-mass radius, R - binary 


ter given in, for example, lHarri 3 d 199(1 updated 2010)) with some 
models having initial concentrations C as large as 125. Roughly 
half of all models had binary fractions fb of 0.1, while fb varied 
from 0.3 to 0.95 in the other half. The adopted binary period dis¬ 
tribution was uniform in the logarithm of the semi-major axis with 
a cut-off at 100AU for all standard models, and either 200AU or 
400AU for the othei^models. For models with large binary frac¬ 
tions, however, the I Kroupa] l i 19951) period distribution was used in¬ 
stead, with a cut-off at log(P / days) = 8.3. Also, we adopted dif¬ 
ferent IMFs for some of ou r models : the Kroupa canonical IMF - a 
two segmented IMF l(Krougajj200l|), the Kroupa stan dard IMF - a 
three segmented IMF ( Kroupa, Tout & Gilmorell993l) . and the two 
segmented modified Kroupa IMF (with different power-law indices 
above and below the break mass). The minimum and maximum 
masses for all IMFs were taken to be O.O8M 0 and 100 .OMq, re¬ 
spectively. A metallicity of Z = 0.001 was adopted for all models. 
Supernovae (SN) natal kick velocities for BHs were modified for 
some mo dels acco rding to the mass fallback procedure described 
bv lBelczvnski. Kalogera & Bulikl 120021) . 


At the end of this Section, we define the terminology used 
from now on in the paper. A hyperbolic collision refers to a physi¬ 
cal collision between two single stars, assuming the sticky-star ap¬ 
proximation (i.e. if the radii of two stars overlap, they are assumed 
to collide and form a single new star). A collision refers to a physi¬ 
cal collision (sticky-star approximation) between two or more stars 
during dynamical interactions between a binary star and either an¬ 
other single star or binary. A merger refers to binary coalescence 
driven by stellar or binary evolution. Mergers can be forced by ear¬ 
lier binary-mediated dynamical interactions. These serve to make 
binaries harder and/or more eccentric, or drive direct collisions. 
Thus, collisions and mergers both involve binaries, but hyperbolic 
collisions involve only single stars. 

The time axis of most of the figures presented in this paper 
extend beyond a Hubble Time - in some models up to 80 Gyr! 
Such exceedingly long evolutionary time scales are irrelevant for 
direct comparison to Galactic GCs. However, they offer important 
insight regarding what to expect for the spatial structure of dissolv¬ 
ing clusters containi ng an IMBH, which could h elp to identify them 
observationally (e.g. lBaneriee & Kroupall201 it) . 


4 RESULTS 

While anal ysing the simulation data for the first project 
dLeigh et alj|2013bl) . it was unexpectedly noticed that a substantial 
buildup of BH mass occurs in some models (see Fig. |6j. We find 
two different regimes of BH mass growth. The first begins at later 
times in the cluster evolution (usually close to the core collapse 
time) and has a small rate of mass accretion (SLOW scenario). The 
second begins near the very beginning of our simulations (usually 
just after the formation of NSs and BHs via supernovae explosions 
and formation of a dense self-gravitating BH sub-system), with a 
very high rate of the mass buildup (FAST scenario). A detailed 
dliscussion of these phenomena is postponed to the next subsec¬ 
tion. We begin by concentrating on the dynamical mechanisms that 
lead to the formation of NSs and stellar mass BHs, which ultimately 
act as the seeds for IMBH formation in the SLOW scenario. 


4.1 Neutron Stars and Stellar Mass Black Holes 

Generally, it is widely assumed that most NSs and BHs in GCs, 
if they are present, formed due to ordinary stellar evolution. This 
leads to SN explosions very early on in the cluster evolution, which 
impart significant natal kicks to NSs and BHs, with velocities rang¬ 
ing from a few to several hundreds km/s. NSs and BHs that receive 
kicks larger than the escape velocity are quickly removed from the 
system. Others remain in the cluster, and experience dynamical in¬ 
teractions with other stars and bin aries that strongly shape the clus¬ 
ter evolution JContenta, Varri & Heggiell2015l) . It is believed that 
new NSs and BHs can form via c ollisions between stars only in 
extremely dense environments (e.g. lPapish & Peretsl2015 . and ref¬ 
erence therein), where the collision rate is sufficiently large. In this 
section, we consider the possibility of forming NSs and BHs in 
less dense environments, not only due to physical collisions during 
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Wo = 6, T. Underfilling, RG=8kpc, M=500000 Mo, bf=0.95, Rh/Rt=50, Rt=98 pc 



Figure 1. The numbers of different kinds of objects containing NSs and/or 
BHs as a function of time, for models in which an IMBH forms. The model 
parameters are given in the figure title. NS - solid line, BH - long dashed 
line, NS-ANY - short dashed line, BH-ANY - dotted line. Here, ANY refers 
to the MS and to all stars that have evolved away from the main sequence. 

single-single, single-binary and binary-binary encounters, but also 
due to the mergers of binaries via normal stellar/binary evolution 
or dynamically-forced channels. 

In Figs. |T| and [2] we show the evolution in the numbers of dif¬ 
ferent kinds of objects containing NSs and/or BHs, both for models 
featuring IMBH formation and those that do not. Evidently, most 
NSs and BHs are quickly evacuated from the system because of the 
large natal kicks they receive (the presented models have a natal 
kick velocity dispersion equal to 265 km/s for both NSs and BHs). 
Only about 20-30 NSs and a few BHs are left, which then begin to 
interact dynamically with other stars and binaries. In both figures, 
the numbers of NSs experiences a sharp initial decrease, followed 
by a slow rise at later times. This is because new NSs are formed 
due to collisions during dynamical interactions and mergers due to 
binary evolution-mediated mass transfer (see Figs. [3] and 0. with 
the latter dominating substantially over the former. 

When the mass of BH is greater than about IOOA/q and IMBH 
mass growth begins (see Fig.Q]- time around of 15 Gyr), all other 
BHs are removed from the system due to strong dynamical interac¬ 
tions with the IMBH. Only one BH remains in the system, slowly 
increasing its mass. Similarly, the number of NSs starts to decrease, 
albeit at a slower rate. If an IMBH does not form, on the other hand, 
a very different situation arises. In these models, the number of NSs 
and BHs (both single and in binaries) is much larger at later times. 
Consequently, a much wider variety of exotic objects containing 
NSs and BHs are present. This leads to an interesting constraint, 
which can be verified observationally - that the number of BHs (or, 
more generally, X-ray binaries) observed in GCs can serve as an in¬ 
dicator of the presence of an IMBH. Systems with an IMBH should 
not contain any stel lar mass BHs. This conclusion is in agreement 
with the findings of iLeigh et al.l d2014bl) . It is worth stressing that 
not all of the NSs and BHs present in our simulations formed via 
SN explosions near the very beginning of the cluster lifetime . They 
also form later on in the cluster evolution, either via collisions also 
hyperbolic collisions) or binary mergers/mass-transfers. 

The radial positions at which NSs and BHs are formed within 
the cluster during its evolution are show in Figs.[3]and0for clusters 
with and without IMBHs, respectively. The BHs are almost exclu¬ 
sively formed at the beginning of the cluster evolution due to stellar 
and binary evolution of massive MS stars. In contrast, NSs are also 


Wo = 6, T. Underfilling, RG=8kpc, M=500000 Mo, bf=0.95, Rh/Rt=50, Rt=98 pc 



Figure 2. The numbers of different kinds of objects containing NSs and/or 
BHs as a function of time, for models in which an IMBH does not form. 
The model parameters are given in the figure title. The various types of 
interactions are described in the caption for Figure[T| 


formed later on in the cluster evolution, but here they are mainly 
formed via collisions during binary-binary dynamical interactions 
or via single and binary star evolution. The NS formation rate is 
similar in both sets of models (i.e. with and without an IMBH) up 
to the time of core collapse or IMBH formation. The times at which 
these NS form, which stretch over many Gyr, are much larger than 
the evolution time for the massive MS stars formed at time T = 
0. Even at these late times, stellar evolution can be rejuvenated by 
physical collisions, which increase the stellar mass. In more com¬ 
plicated ways, binary evolution can have a similar effect: earlier dy¬ 
namical interactions of binaries can make a binary harder, or lead 
to collisions during interactions. Thus these figures implicitly indi¬ 
cate very clearly the complicated interplay between stellar dynam¬ 
ics and stellar evolution. For models in which an IMBH is formed, 
the late-time formation of NSs and/or BHs is strongly inhibited. 
The core radius remains large, and the probability of strong dy¬ 
namical interactions is low. Furthermore, most of the massive ob¬ 
jects (stars or binaries) that act as the progenitors for NSs and/or 
BHs are efficiently removed from the system due to dynamical in¬ 
teractions with the IMBH. A different situation arises in models 
without IMBH formation. Here, the cluster reaches the gravother- 
mal oscillation phase, such that the core bounces in and out of a 
very dense state. This environment is ideal for frequent and strong 
dynamical interactions. Indeed, a substantial number of NSs and a 
few BHs are typically formed during this phase of cluster evolu¬ 
tion, only some of which form via normal stellar/binary evolution. 
Again, early physical collisions and/or binary dynamical interac¬ 
tions stimulate or induce accelerated binary and/or single star evo¬ 
lution. In the case of binaries this can occur via the hardening of 
binaries due to dynamical interactions, provided collisions do not 
happen first. Binary hardening ultimately increases the probability 
of NS/BH formation due to single/binary evolution. Interestingly, 
a significant number of NSs formed via binary evolution are found 
well outside the core, and even beyond the half-mass radius. The 
progenitor binaries experienced a strong dynamical interaction with 
the IMBH deep in the cluster core and were subsequently kicked to 
large clustercentric radii before (or after) NS formation due to sin¬ 
gle/binary evolution occurred. This suggests an interesting obser¬ 
vational constraint for the presence of IMBHs in globular clusters: 
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Time (Myr) 

Figure 3. The position with respect to the centre of the cluster at which NSs 
and BHs are formed for models with IMBH formation. The model parame¬ 
ters are given in the figure title. R c and Rf, are the core and half-mass radii, 
respectively. BIN-SIN: binary-single dynamical interaction during which 
a collision occurs, BIN-BIN: binary - binary dynamical interaction during 
which a collision occurs, COLL: single - single collision in hyperbolic en¬ 
counter, SIN-EVO: single star evolution and BIN-EVO: binary star evolu¬ 
tion. 


Wo = 6, T. Underfilling, RG=8kpc, M=500000 Mo, bf=0.95, Rh/Rt=50, Rt=98 pc 



Figure 4. The position with respect to the centre of the cluster at which 
NSs and BHs are formed for models without IMBH formation. The model 
parameters are given in the figure title. R c and /iy, are the core and half¬ 
mass radii, respectively. The various types of interaction are described in 
the caption for Figure [3 


Clusters observed to harbour a significant fraction of their NSs be¬ 
yond Rh could be more likely to host an IMBH. 

The conclusions described in the above paragraph regarding 
the importance for NS/BH formation via collisions, mergers or bi¬ 
nary mass transfer (due to both dynamical interactions and nor¬ 
mal stellar evolution) are further supported by Fig. [5] These figure 
shows the types of stars which will form NSs or BHs during either 
stellar/binary evolution or dynamical interactions. The various sta r 
types are the same as used in the BSE code l lHurlev, Pols, & Tout! 
l200d) . namely: 0 and 1 - MS, 2 - 9 - evolved stars, 10 - 12 — 
white dwarfs (WDs), 13 - NS, 14 - BH. NSs created in the course 
of the cluster evolution, long after the initial SN explosions of the 
most massive stars, are primarily formed due to binary evolution 
in binaries containing a MS star and an evolved star, containing 
two evolved stars or containing a WD and an evolved star. There 


Wo = 6, T. Underfilling, RG=8kpc, M=500000 Mo, bf=0.95, Rh/Rt=50, Rt=98 pc 
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Figure 5. The different types of interacting stars that will form a NS or BH 
via a merger/collision for models without IMBH formation. The model pa¬ 
rameters are given in the figure title. Initial Type refers to the initial type of 
a star, which will ultimately merge/collide to form a NS or BH. Companion 
Type refers to the type of star acting as a binary companion to the collid¬ 
ing star. The different types of interactions are described in the caption for 
Figure [3 

are also merger events connected with gravitational wave radiation 
(GR) in WD-WD and NS-WD binaries in the case of BH forma¬ 
tion. In order for orbit shrinkage due to GR to be efficient, the bi¬ 
nary semi-major axis must be on the order of a solar radius. Such 
small semi-major axes may result either from pure binary evolu¬ 
tion involving a common envelope stage (about 30% of all WD- 
WD GR mergers), or from binary evolution forced by dynamical 
interactions with other stars and/or binaries (frequently connected 
with the exchange of a MS star into a binary and then followed 
by a common envelope stage). Dynamical interactions can lead to 
gradual binary hardening and substantial eccentricity pumping, un¬ 
til the semi-major axis is sufficiently small and the eccentricity suf¬ 
ficiently large for GR to become efficient. Usually, these dynamical 
interactions are not sufficiently strong to make binaries escape from 
the cluster, but they are strong enough to eject them from the core. 
Interestingly, there is a substantial number of physical hyperbolic 
collisions between a WD and an evolved star and even between two 
WDs. Generally, NSs and BHs that form due to either collisions 
during dynamical interactions or mergers and binary mass transfer 
(connected with stellar evolution) have all types of stars as their 
progenitors. 

4.2 A New Scenario for IMBH Formation 

As already mentioned in Section[4] there are two different regimes 
of BH mass buildup (see Fig. [6] and the ’’white space” - a region 
between about 1.0 to 3.0 Gyr with lowest density of lines, which 
marks approximately the separation between these two regimes). 
One starts later on in the cluster evolution (usually close to the core 
collapse time) with a small rate (SLOW scenario), while the other 
starts almost immediately from the beginning of our simulations 
(usually just after the cessation of NS and BH formation via su¬ 
pernovae explosions and formation of a dense self-gravitating BH 
sub-system) with a very high rate of mass buildup (FAST scenario). 
As discussed in detail later on in this section, the FAST scenario re¬ 
quires that a dense and massive BH subsystem forms shortly after 
the most massive stars explode as SN. The SLOW scenario, on the 
other hand, requires that all BHs formed during SN explosions are 
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Figure 6. IMBH mass buildup for models with IMBH formation. Dashed 
lines - models with reduced mass accretion onto the IMBH (25% of the 
standard BSE setup) and reduced (relative to the default) star expansion af¬ 
ter merger events (1/3 of the standard Fewbody setup). Solid lines - mod¬ 
els with the standard prescription for IMBH mass accretion rates and post¬ 
merger expansion. See description in the text. 


Figure 7. Time evolution of the semi-major axis for binaries containing an 
IMBH in the SLOW scenario. The model parameters are given in the figure 
title. Symbols: plus - flybys or exchanges during binary interactions, cross 
- collisions during binary interactions in which the binary is preserved, star 
-total collisions, where all dynamically interacting stars collide to form a 
single object. 


ejected from the cluster, such that the seed BH forms later due to 
dynamical interactions. A substantial buildup of BH mass is found 
in run s both with a nd witho u t mass fallback during SN explosions 
teelczvnski, Kalogera & Bulikl2002l) . Only a small fraction (about 
20%) of all models considered here show a substantial buildup of 
BH mass and subsequent IMBH formation. In general, the process 
of IMBH formation is highly stochastic. Upon inspecting our simu¬ 
lations, we find the following trends for IMBH formation: the larger 
the initial cluster mass, the higher the probability of IMBH for¬ 
mation; the larger the initial cluster concentration, the higher the 
probability of IMBH formation; the larger the initial cluster con¬ 
centration, the earlier and faster an IMBH will form. 

In what remains of this paper, we address the following ques¬ 
tions: 


(i) What are the mechanisms responsible for IMBH formation? 
- Section l4~2l 

(ii) How is it possible to form an IMBH so late in the cluster 
evolution? - Section l4~2l 

(iii) What are the observational (direct and indirect) signatures 
of the presence of an IMBH in GCs? - Section[5] 


The simulations presented in Fig. [6] (solid lines) were carried 
out with the assumption that 100% of the mass of a star colliding 
with a BH is accreted onto the BH, and after the collision the size of 
the final object is three times larger than the sum of the radii of col¬ 
liding stars (Standard accretion case). The sticky-star approxima¬ 
tion is the standard assumption applied in the McScatter interface to 
the BSE code fHeggie, Porte gies Zwart, & Hurlevlf2006l ) and in the 
Fewbody code l lFregeau et alj2004h . Additionally, the size of a star 
post-collision is set to three times the radius of a normal MS star 
with the mass of the collision product, which is the default assump¬ 
tion in the Fewbody code. The process of mass accretion of any 
type of star onto the BH is very complicated (e.g. jets, strong radi¬ 
ation heating the gas, etc. - all of them not included in simulations 
presented in this paper) and recent simulations suggest that less 
than 50% of the incoming star mass can be directl y acc reted onto 
the B H (see e.g. Guillocho n & Ramirez-Rui^l20 1 3l : Kvutoku et al.l 
120151 : Shiok awa et alj 20151) . Metzger & Stonel d2015l) argue that 
the accreted mass fraction can be as small as 10%, and even as 


small as 1%. Such small fractions for the accreted mass are debat¬ 
able given very short time scales required for SMBH formation in 
the early Universe. Notwithstanding, we stress that if the BH ac¬ 
cretion efficiency is much lower than assumed in this paper, then 
no IMBH should form within a Hubble Time. The exact efficiency 
limit is difficult to quantify given the stochasticity of IMBH for¬ 
mation in our simulations. We caution that the efficiency of IMBH 
formation is sensitive to the maximum initial BH mass, which is 
determined both by the maximum IMF mass and the amount of 
mass lost due to stellar evolution. These issues will be the subject 
of future work. In the next set of simulations, to test the IMBH for¬ 
mation scenario in much less favourable conditions, we therefore 
weaken the standard assumptions - only 25% (and even 10% for a 
few cases) of the mass of the incoming star was accreted and the 
size of the final merged object was just the sum of the radii of the 
colliding stars (Reduced accretion case). As is clear from Fig. [6] 
(see the dashed lines), IMBHs still form. As expected for the Re¬ 
duced accretion case, the IMBH formation process is less efficient, 
but still two main scenarios, SLOW and FAST, are present. Typ¬ 
ically, however, the final IMBH masses are lower in the Reduced 
accretion case. 

Importantly, accreting BHs almost always have a binary com¬ 
panion (whatever its nature). Figs. [7] and [8] show the time evolu¬ 
tion of the semi-major axes of binaries containing an IMBH in the 
SLOW and FAST scenarios, respectively. In the SLOW scenario, 
there are clearly visible binary evolution patterns (see the lines of 
plus signs extending nearly perpendicular to the time axis) - shrink¬ 
age of the binary semi-major axis (binary hardening). This is due 
to dynamical interactions with other binaries and single stars, as 
well as gravitational wave radiation (GR). Binary mergers due to 
GR are rare. More frequent are collisions connected with binary 
dynamical interactions (see Fig.ll9t. In the FAST scenario, binary 
dynamical interactions are primarily connected with collisions or 
total collisions due to the extremely high densities (see Fig. EU), 
small binary semi-major axis (see Fig. [8} and very large IMBH 
masses (see Fig. [6j. The evolution of the binary semi-major axis 
displayed by the short lines of cross signs extending parallel to 
the time axis, visible in Fig. [8] is connected with collisions, not 
with flybys. There are two kinds of dynamical collisions in binary- 
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Wo = 6, N=300000, fb=0.1, rplum=100, Rt=69.0 pc, fallback 



Figure 8. Time evolution of the semi-major axis for binaries containing an 
IMBH in the FAST scenario. The model parameters are given in the figure 
title. Symbols: plus - flybys or exchanges during binary interactions, cross 

- collisions during binary interactions in which the binary is preserved, star 

- total collisions, where all dynamically interacting stars collide to form a 
single object. 


Wo = 6, N=300000, fb=0.1, rplum=100, Rt=69.0 pc, fallback 



Figure 9. Companion masses for binaries containing an IMBH and partici¬ 
pating in dynamical interactions in the FAST scenario. The model parame¬ 
ters are given in the figure title. Symbols: plus - flybys or exchanges during 
binary interactions, cross - collisions during binary interactions in which the 
binary is preserved, star - total collisions, where all dynamically interacting 
stars collide to form a single object. 


single and binary-binary interactions: the collision of an incoming 
star with one binary component (which preserves the binary), and 
the collision of all interacting stars three or four into one object (i.e. 
total collision). During a total collision, there is a sequence of sub¬ 
sequent collisions of a puffy star formed in the first collision with 
other interacting stars. For such an event to occur, at least one star 
must not be an evolved compact star (WD, NS or BH), otherwise 
the collision probability is extremely small. Once again, we stress 
that within the framework of the Fewbody code, any dissipational 
processes are not taken into account. So, there are no GR-mediated 
mergers during dynamical interactions. These can only happen via 
binary evolution, which is taken care by the BSE code. The to¬ 
tal collision events are crucial for the IMBF1 formation. Because 
of these events, an IMBF1 binary (when its mass is still low) can 
avoid hardening to the point of being able to escape from the clus¬ 
ter, which can occur if it receives a strong recoil during a subse¬ 
quent dynamical interaction. So, the IMBFI remains in the system 
and steadily grows in mass. 

In Fig. [9] we show the masses of IMBH binary companions. 
As is clear, in the FAST scenario, the early fast increase of the 
IMBH mass is connected with the formation of a massive BH sub¬ 
system in the central cluster regions, and the subsequent dynami¬ 
cal collisions and total collisions of massive stellar mass BHs and 
massive binaries containing BH.{]. In these models, a large num¬ 
ber of BHs remain in the system due to a reduction in the imparted 
kick velocity that follows from the assumption of significant mass 
fallback onto the BHs. These BHs have masses up to 60 Mq, and 
are usually the companions of binaries with accreting IMBHs. Af¬ 
ter about 2 Gyr, these very massive objects are no longer present 
in the system. They end up being removed by, or colliding with, 
the IMBH due to dynamical interactions. At this time, the IMBH 
companion masses in the FAST scenario are similar to those in 
the SLOW scenario. These companions have masses around 1 Mq, 
which are characteristic of the most massive objects remaining in 


1 The strong recoils due to gravitational radiation upon the mergers of BHs 
is not considered and quantified here. The inclusion of this effect can slow 
down IMBH mass buildup, but not prevent IMBH formation. 


Wo = 6, N=300000, fb=0.1, rplum=50, Rt=69.0 pc 



Figure 10. Companion types for binaries containing an IMBH that partici¬ 
pate in dynamical interactions in the SLOW scenario. The model parameters 
are given in the figure title. Symbols: plus - flybys or exchanges during bi¬ 
nary interactions, cross - collisions during binary interactions in which the 
binary is preserved, star - total collisions, where all dynamically interacting 
stars collide to form a single object. 

the system at this time, such as NSs, WDs, MS stars and even the 
collision products of unevolved stars. So, after all stellar mass BHs 
have been removed from the system, the rate of IMBH mass growth 
substantially decreases, becoming more or less similar to what is 
seen in the SLOW scenario. 

The conclusions reached in the previous paragraph are further 
supported by the data presented in Figs.[T0]and[TT] where the types 
of IMBH binary companion stars are shown. In the SLOW sce¬ 
nario, events involving any type of unevolved companion stars are 
spread out over the whole cluster evolution. In the FAST scenario, 
however, they are present only at the beginning of our simulations, 
with a few additionally events involving asymptotic giant branch 
stars occurring later on in the cluster evolution. In both cases, the 
most frequent type of dynamical interaction involves an IMBH bi¬ 
nary with a WD companion. Interactions with IMBH binaries con¬ 
taining NSs and BHs are in general much more rare, but occur more 














MO CCA Code for Star Cluster Simulations - IV. A new Scenario for IMBH Formation 9 
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Figure 11. Companion types for binaries containing an IMBH that partici¬ 
pate in dynamical interactions in the FAST scenario. The model parameters 
are given in the figure title. Symbols: plus - flybys or exchanges during bi¬ 
nary interactions, cross - collisions during binary interactions in which the 
binary is preserved, star - total collisions, where all dynamically interacting 
stars collide to form a single object. 


Wo = 6, N = 300000, fb = 0.1, Rt = 69pc 

1e+11 

1e+10 

1e+09 

1e+08 
o 

Q. 

| 1e+07 

o 

I 1e+06 
100000 
10000 
1000 

0 5000 10000 15000 20000 25000 30000 35000 40000 45000 50000 

Time (Myr) 

Figure 12. Time evolution of the central cluster density in the SLOW and 
FAST scenarios for IMBH formation. The model parameters are given both 
in the figure title, and within the figure itself. 
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frequently in the FAST scenario (near the beginning of the cluster 
evolution). In the FAST scenario, total collisions occur with the 
greatest frequency, due to the large cluster density (see Fig. ll2t and 
the small IMBH binary semi-major axis (see Fig. 0. The IMBH 
mass becomes so high that gravitational focusing efficiently chan¬ 
nels incoming stars/binaries to drift very close to the IMBH, such 
that strong interactions occur regularly. This serves to trigger addi¬ 
tional merger events with the IMBH. 

As it is clear from Fig. EU the central cluster densities reached 
in the SLOW scenario are not very high - only 10 5 Mq/pc 3 (sim¬ 
ilar to the observed densities of many Galactic GCs). Thus, IMBH 
formation does not require any special sets of conditions. The sit¬ 
uation is different in the FAST scenario. The densities required for 
significant IMBH mass buildup to occur are very high, greater than 
10 8 Mq /pc 3 . These extremely high densities are needed when BHs 
form a bound and very dense subsystem in the cluster centre - col¬ 
lisions of binaries containing BHs must be more efficient than the 
removal of the most massive BHs from the cluster due to strong re¬ 
coils during dynamical interactions. The central escape velocity is 


on the order of a few hundred km/s. Such high densities and central 
escape velocities are unlikely to be reached in the GCs observed 
in the Milky Way0. These densities can be much higher in nuclear 
star clusters (NSCs). Could the FAST scenario for IMBH formation 
discussed in this paper occur commonly in the NSCs of low-mass 
galaxies? 

Based on the results presented in this section, the new scenario 
for IMBH formation proposed in this paper can be summarised as 
follows: 

• There are two possible variants of the IMBH formation sce¬ 
nario: SLOW and FAST. 

- SLOW scenario - either a single BH is left after the early 
phase of SN explosions, or a single BH is formed via mergers or 
collisions during dynamical interactions; 

- FAST scenario - several dozen BHs remain in the system 
after the early phase of SN explosions, and form a dense central 
subsystem. The central density must be extremely high (greater 
than 10 8 Mq/pc 3 ) for an IMBH to form. Alternatively, all BHs 
are quickly and efficiently removed from the system via dynami¬ 
cal interactions. If at least one remains, then the SLOW scenario 
is follower^]. 

• Next, the formation of a BH-binary (i.e. a BH in a binary with 
any type of star as its companion) forms via a three-body interac¬ 
tion. The BH is the most massive object in the cluster, so there is 
a high probability that it will form or it will be exchanged into a 
binary. 

• Dynamical interactions with other binaries and stars: 

- orbit tightening leading to mass transfer from MS/RG/AGB 
companions; 

- exchanges and collisions, leaving the binary in tact; 

- total collisions during dynamical interactions or mergers in¬ 
duced by the emission of gravitational waves - in this case, the 
binary is destroyed and only a BH is left. The single BH then 
forms a new binary via another three-body interaction, which 
is free to undergo subsequent dynamical interactions with other 
single and binary stars, and the process repeats. In this way, the 
BH mass steadily increases. 

It is worth noting that the presented scenarios for IMBH formation 
in GCs, in particular the SLOW scenario, do not require any spe¬ 
cific conditions, unlike other scenarios proposed in the literature. 
IMBH formation occurs solely via binary dynamical interactions 
and mass transfer in binaries. 

Having discussed the two regimes of formation scenarios for 
IMBHs and their conditional requirements, we now turn our atten¬ 
tion to possible observational signatures of the presence of IMBHs 
in star clusters. 


2 T he most dense Galactic GCs have central density around 10®Mg/pc 3 
- see iHarrisl 41996t updated 2010) 

3 Interestingly, a dense, gravitationally bound BH subsystem can form 
without assuming any reduction (due to mass fallback) in th e standard pre¬ 
scription for SN kicks iBelczvnski, Kalogera & Bulikll20021) . A substantial 
number of BHs can be retained in the system post-kick due to the large es¬ 
cape velocity, which can be on the order of about 200 km/s for such dense 
clusters. 
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Figure 13. Surface brightness profiles (for the times given in each inset) in 
the FAST scenario. An IMBH begins forming at the very beginning of the 
cluster evolution. The model parameters are given in the figure title. 


5 OBSERVATIONAL SIGNATURES OF IMBHS IN 
GLOBULAR CLUSTERS 

As already stated in SectionQ] IMBHs are only observable via in¬ 
direct means. What is observed is their influence on the host clus¬ 
ter structure or radio and X-ray emissions connected with accretion 
processes. In the future, it may also be possible to observe GR from 
mergers of compact objects with an IMBH. In this Section, we dis¬ 
cuss possible observational signatures of an IMBH. 

5.1 Surface Brightness and Velocity Dispersion Profiles 

Theoretical investigations and numerical simulations suggest that, 
if an IMBH is present in the centre of a GC, the velocity dispersion 
profile (VDP) should rise toward the cluster centre, and the surface 
brightness profile (SBP) should be shallow with a large core ra¬ 
dius. This suggests that the best candidates to harbour IMBHs are 
massive GCs with large core radii, shallow centralJIBPs and cen¬ 
trally rising VDPs. However, as shown bv lBianchini et all ( l2015h . 
when using integrated light, the observed VDP can be heavily bi¬ 
ased by the presence of a few very bright stars in the cluster core. 
This would obscure any underlying observational signature of an 
IMBH. 

In Figs. [13] and [14] are shown the SBP and VDP for a par¬ 
ticular model featuring IMBH formation. A flat SBP with a large 
core radius and a steeply rising VDP toward the cluster centre are 
clearly visible. The IMBH begins accreting mass almost immedi¬ 
ately, eventually increasing its mass to 7000 Mg. The simulated 
VDPs and SB Ps were used to fit the IMBH mass using Jeans’ 
model (see lLtitzgendorf et alj2013l for a detailed description of the 
method). Fig.[l5]shows that the resulting fit to the MOCCA data is 
quite good. The best agreement is found between about 4 to 12 Gyr 
after t = 0, at which time the IMBH mass is already substantial, 
albeit a small fraction of the total cluster mass. At earlier times, 
the IMBH is quickly growing in mass. At even later times (i.e. > 
12 Gyr), the IMBH mass becomes a substantial fraction of the to¬ 
tal cluster mass (the IMBH mass is about 15% of the total cluster 
mass at t = 20 Gyr). Hence, at early and late times, the applied fit¬ 
ting method reaches its limits: at early times the IMBH mass is too 
small to yield a strong imprint on the observed SBPs and VDPs, 
whereas at late times the IMBH mass dominates the central cluster 
potential. The fact that the MOCCA code is able to correctly repro- 


Wo = 6, N=300000, Rt=69 pc, Rt/Rh=100, fb=0.1, IMBH, fallback 



Figure 14. Velocity dispersion profiles (for the times given in each inset) in 
the FAST scenario. VDP is computed only for stars which My < 6 mag. 
An IMBH begins forming at the very beginning of the cluster evolution. 
The model parameters are given in the figure title. 
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Figure 15. The time evolution of the IMBH mass for a particular MOCCA 
simulation, found from fitting Jeans’ model jLtitzgendorf et a I i 2 01 3ll to the 
simulated SBP and VDP. The model parameters are given in the figure title. 


duce the cluster structure in the vicinity of an IMBH is somewhat 
surprising. This is because MOCCA was not designed to model 
the cluster evolution on time scales relevant to a massive IMBH. 
Hence, the successful application of the MOCCA code to IMBH 
formation highlights a number of strengths inherent to the Monte 
Carlo method for cluster evolution that we now understand. Never¬ 
theless, models that feature IMBHs must be inspected with caution, 
since their presence can be masked by other dynamical effects. 

For example, some models in which an IMBH forms lack a 
centrally rising VDP and/or a shallow central SBP. One such exam¬ 
ple is shown in Figs.|T6]and[T7] In this model, IMBH mass growth 
begins about 10 Gyr into the cluster evolution. Up until this time, 
the SBP shows clear signs of core collapse, namely a centrally ris¬ 
ing profile. Later on, the SBP profile begins to flatten. Meanwhile, 
the VDP does not show the characteristic structure of either a col¬ 
lapsing core or a core hosting an IMBH. Instead, the VDP first 
begins to decrease toward the cluster centre, and then becomes flat. 
Only, at a cluster age of about 50 Gyr is there a clear indication of 
a steeply increasing VDP toward the cluster centre. 

To better understand these effects, we turn our attention to 
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Wo = 6, M=500000 Mo, Rg=4kpc, Rt=62 pc, Rt/Rh=50, fb=0.95, Underfiling 



Figure 16. Surface brightness profiles (for the times given in each inset) in 
the SLOW scenario. An IMBH begins forming at a cluster age of about 10 
Gyr. The model parameters are given in the figure title. 


Wo = 6, M=500000 Mo, Rg=4kpc, Rt=62 pc, Rt/Rh=50, fb=0.95, Underfilling 



Figure 17. Velocity dispersion profiles (for the times given in each inset) 
in the SLOW scenario. VDP is computed only for stars which My < 6 
mag. An IMBH begins forming at a cluster age of about 10 Gyr. The model 
parameters are given in the figure title. 



Time (Myr) 

Figure 18. The distance of the IMBH from the cluster centre along with 
the IMBH mass as a function of time in the SLOW scenario. The model 
parameters are given in the figure title. 


the motion of the IMBH relative to the cluster centre. MOCCA 
treats IMBHs as ordinary objects, albeit extremely massive ones. 
They are not pinned to the cluster centre, as is assumed in Jeans’ 
model. In Fig. |T8] we show the position of an IMBH with re¬ 
spect to the cluster centre as a function of time. Clearly, the mo¬ 
tion of the IMBH is significant, particularly for IMBH masses 
less than a dozen solar masses. There are visible spikes in this 
motion connected with strong dynamical interactions between the 
IMBH and surrounding objects. We now know (e.g. Spurzem, 
Diplom a thesis U niver sity of Goetti ngen 1983, personal communi¬ 
cation, [Zhong^erczi^^S^urzemj yOlJ' and references therein)) 
that motion of IMBHs and SMBHs relative to the centre of mass o f 
their host can destroy a Bahcall-Wolf cusp IB ahead & Wo I (II 19761) . 
which is ultimately responsible for the steeply rising central VDP. 
We stress that within the framework of the MOCCA code, the 
off-centre positions of massive IMBHs can alter the structure of 
the underlying cluster potential, and hence influence the orbits of 
cluster stars. Based on the model presented in Figs. uiinnn 
and an analysis of other models, the presence of an IMBH will 
show a clear signal in the VDP for IMBH masses larger than about 
1000 — 2OOOM 0 . 

5.2 Electromagnetic and Gravitational Radiation from 
Accreting IMBHs 

In Fig. [19] we show the time evolution of the IMBH mass in the 
SLOW scenario for models with Standard and Reduced accretion 
rates (see the definitions of Standard and Reduce accretion in Sec¬ 
tion [4~2t . In the figure caption are given the numbers and proper¬ 
ties of mass-transfer events onto the IMBH, either from a compan¬ 
ion during binary evolution, or because of collisions with incoming 
stars. As expected, the Reduced accretion rate case shows a smaller 
rate of IMBH mass buildup relative to the Standard case. For the 
Standard case, the amount of mass accreted by the IMBH via colli¬ 
sions during dynamical interactions is much larger than for binary 
evolution-mediated mass-transfers or mergers, despite the fact that 
the number of times each event occurs is roughly the same. For the 
Reduced case, the amount of mass accreted by the IMBH is roughly 
the same for both mechanisms, despite the fact that the number of 
binary evolution-mediated mass-transfers or mergers strongly out¬ 
weighs the number of collisions during dynamical interactions. The 
rate of IMBH mass growth due to hyperbolic collisions is much 
larger in the Standard case than in the Reduced case. This is be¬ 
cause the formation of binaries in three-body interactions becomes 
less probable with increasing IMBH mass due to the large acceler¬ 
ations imparted very close to the IMBH. Interestingly, mergers due 
to GR and binary evolution are more frequent in the Reduced case. 
This is connected with the fact that collisions during dynamical in¬ 
teractions produce far fewer total collisions, due to the smaller size 
and hence cross-section of the collided object. Instead, a large num¬ 
ber of GR-mediated mergers are observed, mainly associated with 
IMBH-WD mergers. Because of the large masses of IMBHs, GR- 
induced mergers can be effective for semi-major axes as large as 
several hundred solar radii. In general, collisions, mergers with or 
mass transfer onto an IMBH are potentially observable via the as¬ 
sociated electromagnetic and/or gravitational radiation. However, 
the MOCCA and BSE codes are unable to accurately quantify the 
duration of these events. Obviously, this presents a significant chal¬ 
lenge in estimating the probability of observing them. 

In the FAST scenario (not shown), hyperbolic collisions are 
the dominant mechanism for IMBH mass growth provided the 
IMBH is sufficiently massive (larger than 3000 - 5000 M 0 ). For 
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Wo = 6, N=1800000, Rt=125 pc, Rt/Rh=75, fb = 0.1, no fallback 


Wo = 6, M = 500000 Mo, Rt = 181 pc, Rt/Rh = 50, RG = 10 kpc, fb = 0.95 
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Figure 19. The time evolution of the IMBH mass for models with Stan¬ 
dard (SA) and Reduced (RA) mass accretion rates in the SLOW scenario. 
The model parameters are given in the figure title. Accretion events due 
to binary stellar evolution, collisions and hyperbolic collisions are depicted 
by different symbols: plus signs, crosses and dots, respectively. Also listed 
is the total number of IMBH interactions (Ne). The numbers of particular 
types of interactions, and the total mass gained by the IMBH in each such 
dynamical interaction or binary evolution-mediated mass transfer event, are 
as follows: Mergers with IMBH (RA) - N = 6049, M = 277.1 M©, 
Collisions with IMBH in binary (RA) - N = 1175, M = 236.1 Mq, Hy¬ 
perbolic Collisions (RA) - N = 23, M = 8.1 Mq, Mergers with IMBH 
(SA) - N = 3490, M = 354.1 Mq, Collisions with IMBH in binary (SA) 
- N = 3110, M = 2915.7 Mq, Hyperbolic Collisions (SA) - N = 341, 
M = 254.0 Mq. The numbers of different GR-related events are as fol¬ 
lows: SA case - IMBH-WD mergers = 87, IMBH-NS mergers = 0, IMBH- 
BH mergers = 0, RA case - IMBH-WD mergers = 159, IMBH-NS mergers 
= 12, IMBH-BH mergers - 2. 

smaller IMBH masses, the processes described in the previous 
paragraph apply instead. The large number of hyperbolic collisions 
is connected with the extreme central densities and IMBH masses. 
In the vicinity of a massive IMBH, binary formation during a three- 
body interaction is unlikely, and gravitational focusing increases 
the probability of direct collisions. There are also hundreds of bi¬ 
nary mergers due to GR. 


Figure 20. The escape velocities of objects removed from the system be¬ 
cause of relaxation (star symbol) and dynamical interactions (dot symbol) 
for models that do not feature IMBH formation. The line shows the mass of 
the most massive BH in the system. 
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Figure 21. The escape velocities of objects removed from the system be¬ 
cause of relaxation (star symbol) and dynamical interactions (dot symbol) 
for models featuring IMBH formation - SLOW scenario. The line shows the 
mass of the most massive BH in the system. 
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5.3 Hypervelocity Escapers 


In recent years, a large number of hypervelocity star (HVS) 
cand idates have been reported (just to list the most resent pa- 
pers:|Brown. Geller & Kenvonl2014l : lLi et al.l2012l : |Palladino et al.1 


12014 IZheng et ali 120141 and references therein). More than 50 
stars have now been classified as HVSs with velocities larger 
than 400 km/s. HVSs can obtain their large velocities via vari¬ 
ous channels, including a vari ety of d ynamical mech anisms and 
SN k i cks (e.g. Bau mgardt. Gu alandris & Portegies Zwart 1200a : 
iTaurisI|2015l~) . The dynamical channels that create HVSs are con¬ 
nected with interactions between: a SMBH/IMBH and a binary, a 
SMBH/IMBH-SMBH binary and a (either single or binary) star, 
and a SMBH/IMBH and two unbound stars. The formation of 
IMBHs in our models create HVSs. 

In Figs.[20][2T]and[22]we show the escape velocities of objects 
removed from the system because of relaxation and dynamical in¬ 
teractions, both for models that do not feature IMBH formation 
and for models that do, both in the SLOW and FAST scenarios, re¬ 
spectively. We expect objects that escape because of relaxation to 
have relatively low velocities (kilometres per seconds), and objects 


escaping due to dynamical interactions to have relatively large ve¬ 
locities (dozenths or hundreds kilometres per second). Indeed, this 
is the case for models that do not feature IMBH formation (see Fig. 

However, a more complicated picture emerges for models that 
feature IMBH formation (see Fig.l21landl22t. We observe the fol¬ 
lowing features: when a very massive IMBH forms with a mass 
greater than a few hundred - IOOOMq, it is able to produce relax¬ 
ation escapers with velocities on the order of about 100 km/s. These 
large escape velocities are connected with the very large velocities 
of stars in the immediate vicinity of the IMBH. Clearly, the forma¬ 
tion of IMBHs with masses larger than 100 — 200Mq is connected 
with higher velocities of dynamical escapers. For IMBH masses 
up to 1000 — 2000Mq, a substantial number of escapers are con¬ 
nected with dynamical interactions involving binaries, and produce 
escapers with velocities up to 1000 - 2000 km/s. For larger IMBH 
masses, the number of dynamical escapers decreases, but their ve¬ 
locities increase substantially, up to 10000 km/s or more. This de¬ 
crease in the number of dynamical escapers is connected with a 
decrease in the probability of binary formation. New binaries con¬ 
taining very massive IMBHs rarely form (due to total collisions of 
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Figure 22. The escape velocities of objects removed from the system be¬ 
cause of relaxation (dot symbol) and dynamical interactions (dot symbol) 
for models featuring IMBH formation - FAST scenario. The line shows the 
mass of the most massive BH in the system. 

binaries). Consequently, strong binary dynamical interactions are 
rare, and hence rarely produce escapers. 

Very early on in the cluster evolution, there are also HVS es¬ 
capers (NSs and BHs) connected with SN kicks. 

Concluding, the presence of IMBHs in GCs is potentially in¬ 
directly observable via the presence of anomalously high numbers 
of HVSs with measured 3D space velocities extending back toward 
their former host GCs, which might be identifiable using the GAIA 
satellite. Some of these should be low-mass MS stars (with masses 
close to the turn off mass). 


6 GALACTIC GLOBULAR CLUSTERS 

As orig inally presented in iLeigh et ail d2013bl ) and iLeigh etaf! 
d201 5al) . the GC models discussed in this paper are able to repro¬ 
duce the observed correlations between a number of global GC pa¬ 
rameters, including the total cluster mass, concentration parameter, 
mass function slope and binary fraction. Here, we re-visit this exer¬ 
cise in order to see if models that feature IMBH formation are able 
to produce present-day cluster properties consistent with those of 
real Galactic GCs (Harris 1996, updated 2010). 

We assume that the probability of IMBH formation depends 
mainly on the average binary interaction probability and the es¬ 
cape velocity. The larger the escape velocity, the larger the number 
of NSs and BHs retained in the system, which opens the possibil¬ 
ity for IMBH formation. The larger the probability for binary in¬ 
teractions, the larger the probability of total collisions of binaries, 
which prevent massive BHs from being removed from the system 
via strong dynamical interactions. Hence, we ask: Do models that 
feature IMBH formation occupy regions in Escape Velocity - Prob¬ 
ability of Interaction space different from models that do not feature 
IMBH formation? Naively, we expect models that form an IMBH 
to have a high probability for binary interactions but a low prob¬ 
ability of escape. Conversely, models that do not form an IMBH 
should have a high escape probability but a low probability for bi¬ 
nary interactions. 

The probability of interaction can be approximated by con¬ 
sidering an interaction between an average binary at the soft/hard 
boundary and an average single star in the average cluster envi¬ 
ronment inside Rh . So, for binaries with binding energies equal to 


Figure 23. The interaction probability vs escape velocity at time T=0 for 
all models. Only models for which BHs form with masses greater than 150 
Mq are considered as models containing IMBHs. The filled symbols corre¬ 
spond to models featuring IMBH formation, while the open symbols corre¬ 
spond to models without IMBH formation. The circles correspond to mod¬ 
els that assume the standard accretion rate, while the squares correspond 
to models that assume the reduced accretion rate (see text for more details). 
The line (drawn by eye) is a statistical boundary separating models with and 
without IMBH formation. 

the average single star kinetic energy, the binary semi-major axis 
is a s -h oc Rh/N, where N is the number of stars in the system 
and Rh is the half-mass radius. The average stellar mass is < m > 
= M/N, where M is the total cluster mass. The average cluster 
density and velocity dispersion inside Rh are < rih > oc N/Rh 
and < Vh > oc sqrt(M / Rh), respectively, where < rih > is 
the average number density inside Rh and < Vh > is the average 
system velocity. 

M 1/2 

P s _ h = A——AT (1) 

rTn 

where A is a constant factor and AT is assumed to be 1 Myr. Fig. 
1231 shows all simulated models. Indeed, at T = 0 there is a statis¬ 
tical boundary above which the probability for IMBH formation 
is high (filled symbols). We caution the reader that IMBH forma¬ 
tion is a highly stochastic process. Different statistical realisations 
of models with identical initial conditions can yield very different 
evolutions with or without IMBH formation. The boundary line in 
Fig.[23]is drawn by eye. Models with Standard and Reduced mass 
accretion only partially overlap. For the Reduced case, an IMBH is 
likely to form in models with escape velocities upwards of about 
30-40 km/s. For the Standard case, an IMBH can form in models 
with escape velocities upwards of only 10-20 km/s We redo Fig. 
1231 at time T=12 Gyr in Fig. [24] (using the present-day observed 
and modelled cluster parameters). Here, the picture is less clear. 
Models with and without IMBH formation overlap. We note some 
Galactic GCs occupy regions in Escape Velocity - Probability of 
Interaction space for which MOCCA will produce an IMBH (see 
Fig-1131 1 - These include Omega Cen, 47Tuc, M22 and NGC6293. 

Concluding, it is worth mentioning that models with reduced 
natal kicks (because of mass fallback) for BHs may contain a sub¬ 
stantial number of stellar mass BHs even after a Hubble time. The 
number of retained BHs depends on the total cluster mass and 
conc entratio n, an d ther efore the cluster half-mass relaxation time 
jBreen & HeggiellizOlSl) . The longer the half-mass relaxation time, 
the larger the number of retained BHs. This means that stellar mass 
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Figure 24. The interaction probability vs escape velocity at time T=12 for 
all models. Only IMBHs with masses greater than 150 Mq are considered. 
The filled symbols correspond to models featuring IMBH formation, while 
the open symbols correspond to models without IMBH formation. The cir¬ 
cles correspond to models that assume the standard accretion rate, while the 
squares correspond to models that assume the reduced accretion rate (see 
text for more details). The triangles c orrespond to the observed GC param¬ 
eters taken from the Harris catalogue ( lHarri si 19961 updated 2010) 

BHs are most likely to be observed in massive and low concentra¬ 
tion (large) GCs, instead of massive and dense GCs. Also, simula¬ 
tions of dense and massive GCs show that NSs and BHs can form 
in substantial numbers because of hyperbolic collisions and binary 
collisions during binary dynamical interactions later on in the clus¬ 
ter evolution, not only because of supernovae explosions early on 
in the cluster evolution. Interestingly, IMBH formation suppresses 
the production of binaries with NS or BH companions, and even 
the formation of NSs and BHs themselves. This is because these 
objects (WDs, NSs and BHs) tend to be quickly removed from the 
system via dynamical interactions with the IMBH. These are the 
most massive objects in the cluster, so that the probability of their 
being involved in dynamical interactions is the greatest. 


7 CONCLUSIONS 

The MOCCA code (Giersz et al 2013) is at present one of the most 
advanced numerical codes for stellar dynamical simulations, capa¬ 
ble of following the evolution of real star clusters in detail compa¬ 
rable to N-body simulations, but orders of magnitude faster (about 
a day for IV = 2 x 10 6 ). The ability of the code to model the 
evolution of specific star clusters (M4, NGC6397, 47Tuc and M22) 
and to reconstruct the observed correlations of various cluster pa¬ 
rameters have been demonstrated in previous papers. The new sce¬ 
nario for IMBH formation in GC s discussed in this paper ca me as a 
byproduct of two earlier projects iLeigh et al.l2013blJ2015al) which 
aimed to explain the observed correlations between cluster mass 
and binary fraction and between cluster concentration and mass 
function slope. This new scenario does not require any specific con¬ 
ditions. IMBHs form solely via binary dynamical interactions and 
mass transfer in binaries, with the latter playing an especially im¬ 
portant role by inducing collisions. There are two kinds of dynam¬ 
ical collisions: collisions involving one binary component (which 
preserve the binary), and total collisions involving all interacting 
stars which collide to form a single object. Total binary collisions 
are especially crucial for IMBH formation. This is because IMBH 


binaries (when their masses are still low) can avoid hardening to 
the point of being able to escape from the cluster. This occurs if the 
IMBH binary receives a strong recoil during a dynamical interac¬ 
tion. Instead, the IMBH remains in the system and, consequently, 
is able to steadily grow in mass. 

It is worth noting that to observe growth of IMBH mass one 
needs clusters which consists of hundreds of thousands stars and 
are fairly concentrated. Such conditions are difficult (and for a 
while still infeasible) for direct N-body simulations and therefore 
the process of IMBH formation was not observed yet in N-body 
simulations. The one exception are N-body simulations of dense 
open clusters done by Olczak (private communication). He ob¬ 
served formation of BH with masses of about 200 Mq during the 
first a few hundred Myr of cluster evolution. Such BH mass growth 
was also seen in MOCCA simulations for the same initial model. 
Also, Monte Carle models of star clusters simulated so far were not 
enough initially concentrated and massive. 

We find two different regimes of BH mass growth - the SLOW 
and FAST scenarios. The SLOW scenario is initiated at later times 
in the cluster evolution, has a small accretion rate and requires mod¬ 
est cluster densities of around 10 3 Mq /pc 3 . The FAST scenario is 
initiated practically from the very beginning of our simulations, has 
a very high accretion rate and requires extreme densities of around 
10 8 Mq/pc 3 . Generally, only a fraction (about 20%) of all models 
considered here show significant BH mass growth and IMBH for¬ 
mation. Thus, the process of IMBH formation is highly stochastic. 
The larger the initial cluster concentration, the larger the probability 
of IMBH formation and the earlier and faster the IMBH is formed. 

IMBHs can only be observed indirectly via their influence on 
the host cluster structure, populations of particular types of objects 
(including binaries, WDs, NSs, BHs or hypervelocity star escapers) 
or radio. X-ray and gravitational wave emissions connected with 
accretion processes and/or the mergers of compact objects. 

Below we provide our main conclusions, which summarise the 
paper. 

• NSs and BHs can form (in substantial numbers) in the course 
of star cluster evolution due to dynamical interactions (collisions 
and binary interactions). The number of NSs/BHs (single or in bi¬ 
naries) retained in the system is much larger in systems without 
an IMBH. If an IMBH is present, all NSs and BHs are quickly re¬ 
moved from the cluster via dynamical interaction with the IMBH; 

• If the system density is on the order of 10 5 Mq /pc 3 and only 
one stellar mass BH is present, it is possible that BH mass growth 
will proceed due solely to dynamical interactions and mass transfer 
from binary companions. In this SLOW scenario, IMBH forma¬ 
tion happens at late times in the cluster evolution, typically during 
the post core-collapse phase of evolution. The SLOW scenario for 
IMBH formation is more probable than the FAST scenario; 

• If the system density is extremely high (greater than about 
10 s Mq/ pc 3 ), a gravitationally bound BH subsystem can form 
early on in the cluster evolution. Subsequently, an IMBH can 
rapidly form via dynamical interactions between single and binary 
BHs. In this FAST scenario for IMBH formation, the large cluster 
density (and escape velocity) serves to drive the growth of the most 
massive BH in the system. The FAST scenario might occur more 
commonly in the NSCs of low-mass galaxies than in GCs; 

• The process of BH mass growth and finally IMBH formation 
is highly stochastic. The rate of IMBH mass growth depends sen¬ 
sitively on the cluster density. The larger the density, the higher the 
rate; 

• There are frequent but brief episodes of mass transfer from 
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binary companions onto an IMBH, subsequent mergers and colli¬ 
sions with incoming stars. Therefore, it could be possible to observe 
such events in action via X-ray and/or GR emissions. However, our 
results suggest that the probability of such an observation is very 
low and difficult to estimate given the uncertain nature of the initial 
conditions and the stochastic nature of the evolution; 

• We observe a centrally rising VDP and shallow SBP profiles 
with large core radii in a number of simulations with IMBH for¬ 
mation. For these models, we ha ve checke d that the technique for 
IMBH mass extraction applied in lLiitzgendorf et all d2013l) recov¬ 
ers reasonably well the IMBH masses in our simulations. However, 
a significant fraction of models with IMBH formation yield VDP 
and SBP profiles for which Jeans’ model fitting cannot accurately 
recover the IMBH mass. This is connected with the motion of the 
IMBH relative to the cluster centre, which ultimately erases any 
underlying Bahcall-Wolf cusp and centrally rising VDP profile; 

• Very massive (greater than a few IOOOMq ) IMBHs with bi¬ 
nary companions can produce hypervelocity escapers during dy¬ 
namical interactions, with velocities greater than a few 1000 km/s. 
Most of these escapers are NSs and WDs, but some are MS stars 
with masses near that of the turn off. Observations of these hyper¬ 
velocity stars could provide unambiguous evidence of the presence 
of an IMBH; 

• The present-day properties of models featuring IMBH forma¬ 
tion are consistent with the observed properties of some Galactic 
GCs, e.g. Omega Cen, 47 Tuc, M22, NGC6293. 

The results presented above suggest a new avenue for IMBH 
formation in GCs. Notwithstanding, future studies aiming to repro¬ 
duce these results should bear in mind the following weaknesses 
inherent to the MOCCA code: 

• In the framework of the MOCCA code any dissipational pro¬ 
cesses are not taken into account. The Fewbody code (an important 
ingredient of the MOCCA code) models only direct gravitational 
interactions between stars and binaries. The hydrodynamics, grav¬ 
itational radiation (GR), tidal effects (distortions and dissipation), 
mass-loss, etc. are not accounted for. In general, however, we ex¬ 
pect the inclusion of dissipational effects to increase the rate of 
IMBH mass buildup. 

• In the FAST scenario, we caution that the final IMBH mass is 
sensitive to the assumed initial maximum BH mass, and hence the 
maximum mass, rrimax, of the stellar IMF. We assume m ma x = 
IOOMq for the models presented in this paper, but have checked 
that IMBHs can still form via the mechanism presented here for 
maximum stellar IMF masses as low a few tens of Mg. This caveat 
does not apply to the SLOW scenario, since here the initial BH 
mass is always small, due to forming later on in the cluster evolu¬ 
tion from binary evolution or collisions involving NSs and/or WDs. 

• MOCCA does not cope well with physical processes that have 
a time scale comparable to the local dynamical time scale. Hence, 
for example, the formation and evolution of an IMBH loss cone 
can not be modelled. Having said that, MOCCA does successfully 
reproduce the expected IMBH signatures in the simulated VDP and 
SBP profiles; 

• MOCCA was not designed to follow the dynamical evolution 
of extremely massive objects (larger than a few hundred Mg). This 
is due to the adopted method for computing the cluster potential. 
Nevertheless, the initial IMBH mass growth is modelled correctly 
by MOCCA; 

• MOCCA was not designed to model the relaxation of a clus¬ 
ter hosting a massive IMBH. This is once again connected to the 
adopted method of calculating the cluster potential; 


• The Fewbody code (Fregeau et al. 2004) works well for ex¬ 

treme mass ratios. However, doubts linger regarding Fewbody’s 
prescription for mergers, mass accretion rates and general relativis¬ 
tic effects involving massive IMBHs. Again, we do not expect these 
issues to significantly affect our conclusions regarding the initial 
phase of IMBH growth;_ 

• The accuracy o f the BSE code l lHurlev. Pols. & Tout] l200d ; 
iHurlev. Tout. & Polsl [2002) for modelling binary evolution and 
mass transfer in binaries with extreme mass ratios (and/or ex¬ 
tremely massive compact objects) is suspect. We note, however, 
that the mass transferred onto an IMBH because of binary/stellar 
evolution is not the dominant mechanism for mass growth. Thus, 
IMBH formation occurs even if the contribution from binary mass 
transfer is ignored; 

• Extremely large kicks in the case of mergers of BHs with mis¬ 
aligned spins. Such large kicks can spoil or weaken the FAST sce¬ 
nario, but not influence much the SLOW scenario. 

The issues discussed above except the first one can contribute 
to decreasing the probability of IMBH formation, and/or lower the 
rate of BH mass growth, and should be re-visited in future studies 
using more sophisticated simulation techniques. We emphasise that 
the results presented in this paper are very sensitive to the assumed 
rates of accretion onto BHs, which are highly uncertain. For exam¬ 
ple, our prescription for BH growth does not account for the pos¬ 
sibility that radiation produced during the accretion process could 
ionise the surrounding gas and ultimately drive inverse Compton 
scattering. This could substantially reduce our assumed BH accre¬ 
tion rates. This issue should be re-visited when better prescrip¬ 
tions for BH accretion become available. Regardless, the results 
presented in this paper illustrate that, for our assumed initial con¬ 
ditions, stellar mass BHs that remain in clusters until the present- 
day should be actively accreting for most of the cluster lifetime, 
independent of the actual accretion rates. We have shown that BH 
growth could be significant, and even result in IMBH formation. 
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